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DETERMINATION  OF  THE  SPIN  OF  THE  PROJECTILE 


Abstract 

The  motion  of  the  projectile  during  the  first  period  of  the  yaw  is 
studied,  to  determine  a  criterion  for  the  spin  required  for  a  real  rate 
of  precession.  M.  d'Adhdmar  first  derives  the  general  criterion*: 

S5  +  4  B2  (5"  +  0  )  +4  ABQ^A  =  0.  (30) 

He  rofers  to  M.  de  Sparra's  criteria,  derived  in  different  ways  for  the 
case  of  "perfectly  regular  departure'1 ,  viz.: 


(  >!2)2  >  4  *(Bk  l  )Q,  (3D 

(  i£!)3  >  p(Hl,  1  )„.  (32) 

Ihe  departure  is  defined  as  "perfectly  regular"  if 


6o 


=  f  0 


-  Po  =  Qo  =  r0  =  0. 


Normally,  the  yaw  is  slightly  different  from  zero  and  its  initial 
derivative  6*  is  slightly  different  from  the  regular  value  lTof* 

The  departure  is  then  defined  as  "almost  perfect"  if 


6P  =  5°  t 
2 

*o  =  0, 

Po  "  6o  +  0?> 
o 

q0  -  sin  60, 


(62) 

(63) 


0  See  the  List  of  Symbols  below. 


=  C08  bot 


(64) 

(65) 


r* 

p0  *  S2  +  rw. 

In  this  case,  the  criterion  (30)  takes  the  form: 

S  >  4  B  s(Ekl)  ol  (SO) 

in  which  s  may  "be  greater  than  1/3,  which  is  M.  de  Sparre*s  value  for  the 
case  of  “perfectly  regular“  departure,  as  indicated  by  equation  (32). 

/! 

H.  P.  Hitchcock 
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LIST  OF  SYMBOLS 


i.J 


a 

a 

b 

c 

fl 

fs 

g 

3 

k 

L 

m 

n 

P 

q 

r 

s 

t 

u 

V 

Ax 

Az 

A(t) 

B(t) 

C 


caliber,  or  diameter  of  the  projectile. 

trace  of  GA  on  the  plane  OHN. 

drag  coefficient, 

French  ballistic  coefficient 
(American  C  *  ,00114/c) . 

coP5  . 

uA  6 . 

gravitational  acceleration, 
slope  of  the  rifling, 
coefficient  of  obliquity:  Y/6. 
distance  from  G  to  C. 
mass  of  projectile, 
numerical  coefficient. 

component  of  angular  velocity  of  GUVz  or  Gxyz 
on  Gx. 

component  of  angular  velocity  of  GUVz  or  Gxyz 
on  Gy. 

component  of  angular  velocity  of  Gxyz  on  Gz. 

M.  de  Sparre’s  coefficient. 

time. 

component  of  v  on  OX. 
velocity  of  G. 
f x/(o  (if  constant)  . 
f 2/co  (if  constant), 
p  +  XQ  cos  cot  +  |-k  sin  cot. 
q  +  sin  cot  -  cos  cot. 
center  of  pressure  of  projectile. 
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Fi 

F*. 

F(v) 

G 

GA 

GT 

GAT 

GU 

GV 

GUVz 

GY 

Gxx 

GYi 

Gzx 

GXiYiZ! 

Gx 

Gy 

Gz 

Gxyz 

Hi 

H2 

K 

L 

M 

N 

0 


fi/co  =  P6 . 
f2/oj  =  A  6. 

French  resistance  law. 

center  of  gravity  of  projectile. 

axis  of  projectile,  directed  towards  the  point. 

tangent  of  the  true  trajectory,  directed  in  the 
direction  of  motion  of  G. 

plane  of  yawj  plane  of  resistance. 

axis  1  GT,  related  to  the  projectile. 

axis  1  GTU,  related  to  the  projectile. 

’’related  trihedron”. 

intersection  of  GAT  and  Gxj.yj., 

horizontal  axis,  directed  to  the  left  of  GT. 

axis  in  vertical  plane  containing  GT,  directed 
downwards . 

axis  on  GT,  in  opposite  direction, 
trihedron  with  origin  at  G. 
axis  1  GAT. 
axis  i  GA  in  GAT. 

axis  on  GA,  directed  towards  the  base. 

"free  trihedron”. 

Fi’  (if  constant). 

Fg *  (if  constant). 

radius  of  gyration  of  projectile  about  GA. 

6^/5.. 

ci)^6  +  01^. 
co^  -  <0^6 . 

point  on  GT:  distance  GO  =  1. 
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OH 

ON 

OHNT 

0 

OX 

OY 

OZ 

OXYZ 

P 

Q 

R 

S 

T 

Ui 

U2 

Vo 

Vm 

Wx 

WB 

Zo 

A 

B 


horizontal  axis,  directed  to  the  right  of  GT» 

axis  in  vertical  plane  containing  GT,  directed 
upwards . 

trihedron  with  origin  at  point  0,  and  axes  || 
to  those  of  GXiYxZi  but  oppositely  directed. 

origin  of  trajectory. 

horizontal  axis  in  plane  of  fire,  directed  towards 
the  target. 

vertical  axis,  directed  upwards. 

horizontal  axis  1  OXY,  directed  to  the  right. 

absolute  trihedron. 

Rkz/4  £2  =  M/A  Q  6. 

Jx. 

stability  function:  P/)t»|. 
resistance;  air  pressure. 

( A  8  -  45  Rk l 
tot. 

X0  cos  cot. 

Xc  sin  cot. 

initial  velocity  of  G. 
minimum  velocity  of  G. 

X0  cos  cot  +  sin  cot. 

X  0  sin  cot  -  cos  cot. 

zone  of  trajectory  from  t  =  0  to  t  =  2«/co. 
axial  moment  of  inertia  of  projectile, 
transverse  moment  of  inertia  of  projectile, 
moment  of  external  forces  about  Gx; *Rk  l  6. 
moment  of  external  forces  about  Gy. 
moment  of  external  forces  about  Gz, 


instantaneous  center  of  precession 
retardation:  R/m. 


angle  of  projection:  =  qo. 

angle  between  the  direction  of  resistance  and  Gz. 

yawj^.ZiGz;  tan  6  =  0  o. 

component  of  6  on  OH;  -  6  sin  y.* 

component  of  6  on  ON;  6  cos 

d7 

drift  angle;  tan  r\  =  - 

angle  of  inclination  of  the  tangent  of  the  true 
trajectory,  GT. 

p  cos  cot  +  q  sin  cot, 

p  sin  cot  -  q  cos  cot. 

component  of  angular  velocity  of  GUVz  on  Gz. 

angle  of  inclination  of  the  tangent  of  the  plane 

dY 

trajectory;  tan  t  = 

ratio:  |y'|/|0»|. 

'  ZxGU. 

orientation  of  plane  of  yaw;  precession; 

ZNO  a  ~  ZXiGx. 

A  Q/B . 

component  of  9*  and  n*  on  Gxx. 

component  of  91  and  7}1  on  Gyi . 

component  of  9*  and  V  on  Gzjl. 

component  of  cox  and  co2  on  GY. 
component  of  cox  and  co2  on  Gx. 
component  of  co^  and  on  Gy, 
component  of  co^  and  co^  on  Gz. 


O  0 

A 

$ 

G 

% 


9n  cos  y  +  rjn  cos  9  sin  if  « 

H  y /A  £  6 . 

9n  cosy  +  f)1’  cos  9  sin  Y  -  >1*9’  sin  9  sin  Y 

Spin  of  projectile,  impressed  by  the  rifling. 

p0j  initial  angular  velocity  of  projectile 
about  GA. 
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THE  GYROSCOPIC  MOTION  OF  STABLE  PROJECTILES 

(THIRD  REPORT)* 

DETERMINATION  OF  THE  SPIN  OF  THE  PROJECTILE 

by 

R.  d’Adhemar, 

Engineer  of  Arts  and  Manufactures,  Doctor  of  Sciences. 
(Memorial  de  l’artillerie  fran^aise,  Vol.  12,  1933,  p.  249). 

1.  INTRODUCTION.  Let  Q  be  the  spin  proper,  impressed 
by  the  rifling;  we  have: 

Q  =  2JV0/a; 

Ve  is  the  initial  velocity  on  the  trajectory, 

a  is  the  diameter  of  the  projectile, 

j  is  the  tangent  of  the  angle  of  inclination  of  the 
grooves . 

For  example,  if  a  =  j  =  1/10,  then  Q=  2V0. 

If  the  number  is  not  very  large,  the  pendular  or 
gyroscopic  phenomenon  does  not  exist;  then,  the  drifts  are 
not  regular.  However,  if  the  number  Q  is  too  large,  there 
is  an  inconvenience**:  for,  on  the  descending  branch  of  the 
trajectory,  the  stability  function  Q  may  assume  too  small 
values . 


I  shall  not  study  this  question  in  all  its  ampli¬ 
tude,  and  I  shall  only  seek  a  lower  limit  of  Q. 

M.  de  Sparre  obtained  a  lower  limit,  in  his  Report  *** 
of  1904,  by  using  his  formulas,  which  represent,  approximately, 
the  gyroscopic  motion. 


*  See:  First  report,  Mem.  de  l^art.  fr,  Vol.  VIII,  Ho.  3.  1929 
(Trans.:  Report  No. 72,  Ballistic  Lab. ).• Second  report,  Mem,  de  l'art,  fr. 
Vol.  XI,  No. 3.  1932  (Trans.:  A-IV-42,  Ballistic  Sec.  file). 

**  P.  Charbonnier.  The  gyroscopic  motion  of  the  projectile,  Mem. 
de  1  *  art .  fr. ,  Vol  .VI  . No  .3^1927. 

***  M.  de  Sparre.  The  motion  of  oblong  projectiles  about  their  cen¬ 
ters  of  gravity  ,  Arklv  for  Materaatik.  Astronoml  och  FyslV,  Stockholm,  1904. 
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M.  Esclangon  found  the  true  basis  of  this  discussion*, 
by  starting  from  the  fundamental  equations  of  gyroscopic  motion. 

This  method  offers  perfect  safety.  On  the  contrary, 
there  exists  some  uncertainty  in  the  use  of  an  approximate 
solution,  which  carries  all  the  weight  of  an  accumulation  of 
simplifications;  the  simplifications  may  engender  deformations. 
However,  M.  Esclangon* s  remarkable  analysis  must  be  completed, 
because  he  adopted  the  ballisticians’  first  method  of  approxi¬ 
mations,  as  M.  de  Sparre  did  otherwise  in  his  classic  report  of 
1904«  This  first  approximation,  which  is  always  useful  if  one 
does  not  ask  himself  what  it  gives,  consists  in  the  substitution 
of  the  conventional  yaw  for  the  true  yaw.  Roughly,  it  neglects 
the  curvature  of  the  horizontal  projection  of  the  trajectory; 

I  must  give  a  precise  indication  on  this  subject. 

2.  THE  ROLE  OF  THE  CURVATURE  OF  THE  HORIZONTAL  PROJECTION 
OF  THE  TRAJECTORY.  I  keep  the  notations  of  my  former  reports**: 
GT  is  the  tangent  of  the  true  trajectory,  GA  is  the  axis  of  the 
projectile.  The  plane  OUN  is  perpendicular  to  the  tangent.  The 
axis  OH  is  horizontal  and  directed  towards  the  rear  of  the  plane 
of  the  figure;  the  axis  ON  is  in  the  vertical  plane  containing 
the  tangent  and  directed  upwards.  Let  a  be  the  trace  of  GA  on 
this  plane  of  reference,  so  that: 

0  a  -  tan  5*  6,  GO  =  1 . 

The  coordinates  of  a  are  6X  on  OH  and  6a  on  ON.  I 
designate  by  ip  the  angle  NOa  ,  this  angle  being  counted  counter¬ 
clockwise,  so  that  OH  is  brought  to  ON  by  a  rotation  of  90°. 

Then: 

6X  =  -  6  sin  V  ,  62  -  6  cos  V . ■ 


I  assume  that  the  rifling  is  left  handed:  for  a  well- 
designed  projectile,  the  drift  will  be  to  the  left,  i.e.,  in 
front  of  the  plane  of  the  figure.  This  is  the  plane  containing 
the  plane  trajectory,  which  we  know  how  to  calculate. 

The  drift  angle  is  R  (Fig.  2).  Let  p’  be  the  deriva¬ 
tive  of  R  with  respect  to  time,  t.  The  sign  of  P  arid  that  of  U* 
depend  on  the  sense  of  orientation  that  is  chosen;  but  we  have, 
as  the  absolute  value: 

u  h'l  =  R  (k-l)  I  6X|  .  (1) 


*  E.  Ecclangon.  Motion  of  projectiles  about  their  center  Of  gravity, 
MoiA.  de  1  (art.  fr.,  Vol.  VI,  1927.  No.  3. 

**  R.  dlAdhcmar.  On  the  gyroscopic  motion  of  stable  projectiles, 
Mem,  do  l*axt.  fr.,  1929  and  1932. 


The  resistance  is  designated  by  R  and  the  retardation  by 
R,  so  that: 

R  =  irJ? . 

Replacing  by  its  mean  value  or  else  by  1/Q  (see 
my  previous  reports),  we  have: 


v2  |  nM  =  g  4_Q  k . -  1.  (2) 

m  k  l 

It  is  evident  that  formulas  (1)  and  (2)  are  not  iden¬ 
tical.  For  example,  according  to  (2),  could  not  vanish, 

while,  according  to  (l) ,  h »  vanishes  at  the  origin  if 
does.  This  shows  the  deformation  which  an  element  may  under¬ 
go  as  the  result  of  taking  an  average. 

In  certain  cases,  it  is  necessary  co  make  a  distinction 
between  an  exact  element  x,  and  the_same  element  smoothed 
by  averaging.  The  notations  x  and  x  might  then  be  used.  We 
shall  write,  then: 

^  Vs  In' j  =  g  4A  k  -  1  .  (3) 

.m  k  l 


In  every  question  that  demands  great  precision,  it  is 
the  element  x  which  must  be  studied,  and  not  x.  Here,  for¬ 
mula  (3)  will  be  sufficient. 

Let  us  consider  now  the  angle  of  inclination  of  the  tan¬ 
gent  to  the  trajectory.  For  the  true  trajectory,  this  angle 
is  9.  For  the  plane  trajectory,  this  angle,  at  the  same  in¬ 
stant  t,  is  t  .  If  the  yaw  6  Is  quite  small. 


e  *  x,  '  9»  *  x* . 


I  do  not  actually  insist  on  the  degree  of  approximation  of 
these  formulas,  and  I  recall  that  : 

.  _  fi  cos  T- .  (4) 

dt  ~  r  - v 

so  that  we  can  write:  , 

5«  =  - « ;os- T .  (5) 
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Let  us  form  the  ratio: 

„  =  ULL  =  L  1 21  k  -  1  Ve  . 

I  7=r*  i  m  cos  t  a  k  /  v 


The  ratio  A  /m  is  very  small. 

The  parameter  k  is  poorly  known;  but  we  may  assume  that  its 
value  remains  near  1  or  2  for  very  small  yaws.  This  value  could 
be  3,  4>  5,  ...  for  large  yaws,  in  the  dangerous  zone. 

The  parameter  l  is  the  distance  from  the  center  of  gravity 
to  the  center  of  pressure  C.  If  this  point  C  is  at  the  height  of 
the  middle  of  the  ogive,  then,  in  the  formula  l  =  na,  n  can 
assume  the  values  2,  3,  ... 

Altho  the  ratio  A  /m  is  very  small,  if  cost  is  very  small, 
and  if  the  ratio  VD/v  becomes  quite  large,  it  is  quite  possible 
that  u  is  not  negligible;  the  derivatives  R*  and  Q*  may  have 
the  same  order  of  magnitude. 

In  his  classic  report  of  1904  (p.  237),  M.  Sparre  assumed 
intuitively  that  the  ratio  u  is  always  negligible.  However, 
that  is  not  exact  for  an  almost  horizontal  trajectory  with  a 
very  high  velocity.  V/hat  ballisticians  call  the  "Gamier  ef¬ 
fect”  Is,  we  might  say,  the  fact  that  the  ratio  u  Is  not  negli¬ 
gible,  in  general.  This  remark  has  been  made  previously  by 
Engineer  General  Maurice  Gamier. 

In  the  recent  works  of  M.  de  Sparre  and  Engineer  General 
Charbonnier,  this  ratio  was  no  longer  taken  as  negligible  a  prior 

Prudence  demands  that  we  take  account  of  the  curvature  of 
the  horizontal  projection  of  the  trajectory,  altho  Its  effect 
may  be,  in  certain  circumstances,  very  small. 

3.  FIVE  INSTANTANEOUS  ROTATIONS.  I  take,  as  the  sense  of 
orientation  of  trihedrons,  the  direct  sense  of  trigonometry, 
which  Is  counterclockwise.  G %x  is  the  tangent,  in  the  direction 
opposite  to  that  of  the  motion  of  G.  Gxx  or  GHX  Is  parallel  to 
the  axis  OH,  which  has  been  defined.  Gz  is  the  axis  of  the  pro¬ 
jectile,  towards  the  base;  the  point  would  be  on  the  prolongation 
of  Gz  (Fig.  4). 

The  moment  M  of  the  resistance  defines  the  axis  Gx  (see 
my  previous  reports).  The  axis  Gy  carries  the  second  moment  My 
of  the  forces  of  resistance.  If  it  is  the  Magnus  effect  which 
is  considered,  My  <  0.  If  it  is  the  Esclangon  effect,  My  >  0. 
According  to  M.  Esclangon,  this  couple  My  is  caused  by  lateral 
frictions  and  its  existence  appears  certain  for  a  well  designed 

projectile.  We  have: 
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We  may  write: 


) 


*  Rk  l  6 . 

Jk 


x  -  < 

The  symbol  P  represents  a  known  element.  To  avoid  all 
discussion  on  the  physical  definition  of  M  ,  I  shall  write: 

J 

M  -  ^£2A6  , 

y 

and  I  shall  make  the  apparently  natural  suppositions: 

1°  A/P  is  constant  and  very  small  if  6  is  very  small; 
for  example,  this  ratio  is  equal  to  2%  or  . . . . 

2°  A/P  may  assume  a  more  appreciable  value  if  6  be¬ 
comes  relatively  quite  large;  for  example,  for  yaws  of  about 
12°  or  15°. 

So,  at  the  origin,  we  shppose  A/P  constant  and  very 
small. 


We  have  thus  defined  two  trihedrons,  Gx1y1z1  and  Gxyz. 

Let  GU  and  GV  be  two  axes  in  the  plane  Gxy,  related  to  the 
projectile.  We  have  a  third  trihedron  GUVz,  related  to  the 
projectile.  I  call  it  the  ’’related  trihedron”,  and  the  tri¬ 
hedron  Gxyz  shall  be  the  ”free  trihedron”  (Fig.  4) • 

If  GXjyiZx  were  fixed,  we  would  simply  have  to  consider 
the  three  Eulerian  angles  cp,  v  ,  6  (I  must  substitute  6  for  the 
classic  notation  9,  since  I  have  already  used  the  symbol  G) . 
Furthermore,  we  would  have  only  the  three  instantaneous  rota¬ 
tions: 


dig 

dt 


or  cp’ , 


df 

dt 


or 


vf 


dfi 

dt 


or  6’ 


We  have,  besides,  two  instantaneous  rotations: 

1°  The  rate  of  depression  of  the  tangent,  9’, 

2°  The  angular  velocity  of  drift,  p*  (Figs.  2 
and  3) . 

I  suppose  essentially  that  the  true  trajectory  and  the 
plane  trajectory  are  very  close,  in  size  and  form,  so  that,  at 
each  instant  t; 

9  «  t,  9  ’  »  t  1 ,  9”  *s  t”  . 

I  do  not  actually  discuss  this  approximation,  Y/hich  is 
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based  on  physical  intuition  of  the  phenomenon.  I  shall  always 
designate  the  initial  velocity,  on  the  trajectory,  by  V0,  and 
the  minimum  velocity  by  V^. 

It  must  be  remarked  that: 


In  the  dangerous  zone  of  the  trajectory,  i.e.,  in  the 
neighborhood  of  the  point  where  Q  is  a  minimum  and  in  the 

neighborhood  of  the  point  where  v  is  a  minimum,  the  ratio 

| t  «  [/ | t»| ,  will  be  extremely  small.  On  the  other  hand,  at 
the  origin,  this  ratio  is  not  very  small.  In  fact,  the  ratio 

|vf  |  /v  can  take,  at  the  origin,  valves  such  as  1/10  or  2/10  ... 
So,  at  the  origin,  the  ratio  I'm  |  /  |Tt  |may  be  about  l/lO  or 
2/10  ....,  it  may  not  be  negligible. 

The  sense  of  orientation  is  counterclockwise  (Fig.  2); 
so,  the  axis  0Z  is  directed  towards  the  rear  of  the  plane  of 
the  figure,  the  plane  OXY.  For  an  observer  placed  on  0Z, 

the  angle  T  is  positive  on  the  ascending  branch,  and  negative 

on  the  descending  branch. 

An  observer  placed  on  0Z  sees  the  tangent  turn  in  a  clock¬ 
wise  direction.  So,  this  rotation  is  represented  by  a  vector 
|  0f  |  ,  directed  forwards.  Consequently,  the  component  on  the 
axis  Gxx  is  negative;  it  is  0*,  or  approximately  . 

When  I  supposed  that  the  true  trajectory  is  very  near  the 
plane  trajectory,  I  implicitly  assumed  that  the  projectile  is 
stable,  that  5  remains  small.  Consequently,  the  drift  will  be 
normal;  q  and  q’  will  be  positive  (Figs.  2  and  3). 

Moreover,  the  passage  from  (2)  to  (3)  would  also  imply 
that  the  projectile  is  stable.  If  V0  is  large  and  if  V  ,  the 
minimum  of  v,  is  not  small,  the  drift  will  remain  feeble, 

(k  -  l)  will  not  become  large,  and  q>  will  remain  small  like 


I  shall  not  actually  discuss  the  second  derivative. 


Tin. 


If  the  departure  of  the  projectile  is  quite  regular,  H £  =  0. 
I  designate  an  initial  value  by  x0  or  x°. 
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4.  SOME  INDICATIONS  ON  THE  NUMERICAL  DATA.  Before  going 
farther,  I  shall  give  some  indications,  some  figures,  whose 
knowledge  is  indispensable,  if  we  wish  to  avoid  vague  discussions. 
We  designate  by  A  the  axial  moment  of  inertia;  in  the  formula 
A=  na5  t  n  may  vary,  e.g.,  between  200  and  300*.  We  designate 
by  B  the  transverse  moment  of  inertia;  in  the  formula  0  =  n  A  , 
n  may  vary  from  4  to  8,  approximately. 

In  the  formula  A  -  mK2,  we  see  intuitively  that  K  must 
remain  between  a/2  and  a/3,  nearly. 

If  we  take  J  =  1/10,  that  corresponds  to  an  angle  of 
inclination  of  the  rifling  of  5°43fJ  this  is  a  probable  datum. 

For  an  angle  of  7°,  J  -  0.123.  For  an  angle  of  10°, 

J  =  0.176. 


The  symbol  P,  which  plays  an  interesting  role,  represents 
an  approximate  mean  value  of  the  angular  velocity  of  the  point 
g.  about  the  instantaneous  center  of  precession  p  (see  my  pre¬ 
vious  reports).  We  have; 

P  =  ~ R  xl  ,  and  R  =  cF(v)  . 

Let  us  seek  an  indication  of  the  numerical  values  of  P, 
at  the  origin,  supposing  a  =  1/10  (100  mm  caliber). 

Let  us  take  the  resistance  lav/  in  the  form  bv2. 

If  the  velocity  v  varies  from  400  ro/s  to  1500  m/s,  we  can 
set  b  =  1/3**.  This  is  a  rough  mean,  but  it  is  sufficient  for 
us.  Let  us  take  c  =  1/2000***;  we  are  at  the  origin,  so  we  do 
not  have  to  consider  the  variations  of  altitude.  We  may  suppose, 
at  the  origin,,  that  k  is  about  1  or  2.  Also,  we  can  take  l  =  2a 
or  3a. 


So,  if  we  set 

ki  =  4/10,  K2  =  a2/5; 


•Hera  a  is  the  caliber  in  meters,  and  A  is  the  axial  moment  of  inertia  in 
kg.ra.  divided  by  the  gravitational  acceleration  in  m/sec.2  —  H.P.H. 

♦•Hence,  bv2  is  the  drag,  in  kg.m/soc.2,  of  a  1  m.  projectile,  in  air  with  a 
density  of  1  kg/m^.  This  value  of  b  corresponds  to  a  resistance  coefficient 
Cr  =  12x10~5.  - H.P.Ho 

•••This  value  of  the  French  ballistic  coefficient  c  corresponds  approximately 
to  an  American  ballistic  coefficient  C  -  2.28.  —  H.P.H. 
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we  will  have 


P0  =  Vo/60. 

For  VG  =  600,  P0  =  10. 


If  we  take  c  =  1/4000,  keeping  the  values  of  all  the  other 
parameters,  we  will  have  PG  =  Vo/l20.  For  V0  =  1000,  PG  *8. 

Let  us  also  seek  an  indication  of  the  numerical  values  of 
P*/P  at  the  origin. 

We  suppose  VG  >  400  m/s.  We  can  neglect  the  variations  of 
b,  k,  l  and  take  b  =  1/3,  since  we  are  not  seeking  any  great 
precision  here. 

I  recall  the  general  formula: 

=  vr  =  -g(sin  t  +R  / g)  . 

When  V0  is  a  large  number,  the  ratio  R/g  can  assume  values 
such  as  40  or  50,  e.g.,  or  even  60;  then 


Then,  we  can  write: 

(P»/P)o  *  2(vf/v)o  *  -2bcV 0, 
and  the  approximation  is  good  if  V0  is  large. 

For  c  =  1/2000  and  V0  =  500,  (P»./P)0  =  1/6. 

For  c  =1/4000  and  V0  =  800,  (P»/P)0  =  8/60. 

For  c  =  1/10,000  and  V0  =  1500,  (P«/P)0  =  1/10. 
So: 


(P«/P)o  =  n/10, 

with  a  coefficient  n  about  unity. 

Let  us  return  to  the  formula  for  £2  .  For  the  usual  artil¬ 
lery  calibers,  if  Q=  nV0,  the  coefficient  n  is* between  1  and  4« 
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It  is  often  said  that  S  is  a  large  number.  To  be  precise: 

Q  is  large  with  respect  to  unity,  with  respect  to  P,  and  also 
with  respect  to  |t'|,  if  the  minimum  velocity  V  ,  on  the  trajectory, 
is  not  very  small;  this  exceptional  case  is  excluded. 

In  the  theory  which  is  going  to  be  developed,  the  approxi¬ 
mation  is  very  good  if  the  number 

0)  =  AS  /B  , 

which  is  between  1/6  and  1/8  of  £2  ,  is  still  quite  large  with 
respect  to  P0. 

5.  THE  INSTANTANEOUS  ROTATION  OF  THE  RELATED  TRIHEDRON 
AND  THE  INSTANTANEOUS  ROTATION  OF  THE  FREE  TRIHEDRON.  Let  us 
consider  the  instantaneous  rotations  of  the  two  trihedrons  which 
have  been  defined,  and  let  us  seek  their  components  on  the  axes 
of  the  free  trihedron  Gxyz.  Let  us  first  take  the  related  tri¬ 
hedron.  The  components  of  the  rotational  velocities  6’,  y’ ,  and 
cp’  are: 

6’,  on  Gx;  y  ’  sin  6,  on  Gy;  cos  6  +  ,  on  Gz. 

The  components  of  the  rotational  velocities  9’  and  are: 

(ot  -  9>  <  0,  on  Gxj ; 
to2  =  t}1  cos  9,  on  Gyij 
0)3  =  -  sin  9,  on  Gzx. 

Let  us  consider  the  trace  GY  of  the  plane  of  resistance  on 
the  plane  Gx^x*  The  axes  Gx  and  GY  are  orthogonal,  as  are  Gxx 
and  Gy!  (Fig.  5).  Let  us  take  the  components  of  cox  and  o>2  on 
these  new  axes  Gx  and  GY.  They  are: 

0)3  =  (Ox  cos  Y  +  002  sin'}'  ,  on  Gx; 

o)^  =  -  cox  sin  y  +  o)2  cos  y  ,  on  GY. 

Now,  we  must  take  the  components  of  (o_  and  co,  on  the  axes 
Gz  and  Gy;  this  gives:  4 

cd^  =  0)3  sin  6  +  (o^  cos  6,  on  Gy; 

( Oy  =  0)3  cos  6  -  co^  sin  6,  on  Gz. 
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So,  the  components  of  the  instantaneous  rotation  of  the 
related  trihedron  are: 

on  Gx:  p  =  6’  +  co,., 

or  p  =  s»  +  gt  cos  y  +  q'  ■  cos  9  sin  V; 

on  Gy:  q  -  yf  sin  6  +  oi^, 

or  q  -  y’  sin  6  -  rp  sin  9  sin  6 

+  cos  6  (-9 f  sin  y  +  71*  cos  9  cos  y); 

on  Gz:  p  =  cp1  +  y  f  cos  5  +  a>y, 

or  p  =  q>»  +  y»  cos  6  -  q  f  sin  9  cos  6 

-  sin  6  (-9 1  sin  y  +  q»  cos  9  cos  y  ). 


If  the  yaw  is  quite  small,  we  can  substitute  5  for  sin  6 
and  1  for  cos  6.  Then,  let 


M  =  (0-6+  co. 
i  4 


(7) 


=  -  qT  6  sin  9  +  (-9*  sin  y  +  cos  9  cos  y  ), 


N  =  a>-  -  oo.  6 
->  4 


(8) 


=  -  qf  sin  9-6  (-9*  sin  y  +  q'  cos  9  cos  y  ). 

We  will  have  the  following  expressions: 

p  =  6'  +  9f  cos  y  +  qf  cos  9  sin  y  ,  (9) 

q  **  y*  6  +  M,  (10) 

p  w  (J)l  +  yf  +  N.  (ll) 

To  pass  from  the  related  trihedron  to  the  free  trihedron, 
it  is  sufficient  to  set 


whence 


<P  =  0,  <p»  =  0; 


r  *  y1  +  N. 


(12) 


The  components  of  the  instantaneous  rotation  of  the  free 
l^hedron  areP,  9,  r. 
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I  substitute  M  for  cd^  and  N  for  a>y,  to  simplify  the  ex¬ 
planation;  but  these  approximations  are  not  indispensable. 
Notice  that; 


**  *(,- 

and 

dM  dco/  . 
.«  at  ' 

S  *  eij 

and 

dN  M  do>7 

dt  at2 

In  fact. 

d  s  in  5 

=  cos 

fi  d6  «  d5 

dt 

dt  dt. 

d.  cos  5  =  ~  sin  5  dS  n 
dt:  dt*  °» 


Slnc«£.6 - Tetaaisas  qui te  small,  and  the  smaller  6  is,  the  more  sat 
isfa^tory  .the'  approximation  is.  As  far  as  p  is  concerned,  equn 
tioh ;{?)  is  ex&c t ,  so  that  we  obtain  immediately  by  differr.n 
tial 

Now,  I  an  going  to  define  "perfectly  regular  departure". 

I  shall  say  that  the  departure  is  perfectly  regular  if 

5„  =  V0  *  p»  -  q6  =  To  =  0. 

1  exclude  the  case  of  vertical  fire:  a  =  90°. 

Consequently:  5®  =  0. 

Since  the  initial  yaw  vanishes,  the  gyroscopic  motion  con¬ 
sists  of  oscillations  from  the  origin  to  the  summit,  and  even 
beyond  that  {see  my  previous  reports).  So  v:e  certainly  have, 
at  the  origin,  a* <  0,  and  since  the  sense  of  orientation  is 
counterclockwise,  we  must  write: 

-  u  q*  =  £(k  -  1)5! .  (13) 

On  the  other  hand,  we  would  have  to  take  the  +  sign 
if  the  sense  of  orientation  were  clockwise,  i.e.,  if  GZ  were 
directed  forward  (Fig.  2). 
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('  V 


Since  u0  >0  and  &°  =  0,  whatever  the  value  of  k  may  be, 

at  the  origin,  1 

S"o  =0. 

On  the  other  hand,  with  9  and  t as  defined, 

v9*  =  -  g  cos  9  +R  (k-l)52,  (14) 

with 


V  Tl  *  -  g  COS  T  , 

and,  if  the  departure  is  perfectly  regular,  =  0. 
So,  at  the  origins 

90  =  =  a  and  9£  =  x£, 

while,  in  general, 

9  *  r  and  9f  88  t1  . 

Since 


Po  =  Vo  -  *]£  =  0, 

therefore 

&o  +  To  “  0. 

The  term  T] »  cos  9  sin  y  is  extremely  small  at  the  ori¬ 
gin,  for  it  contains  two  factors  which  vanish.  It  might  be 
said  that  this  term  is  of  the  "second  order";  but  it  is  some¬ 
times  dangerous  to  extend  to  the  "very  small",  notions  which 
are  clear  only  for  the  "infinitesimal". 

Let  us  retain  the  formula; 

A&.+  n  =  or  =  in,  I  (15) 

relative  to  the  case  of  perfectly  regular  departure.  Notice, 
besides,  that,  in  this  cases 


(D»  =  (D»  =  (i)g  =  =  M°  =  N°  =  o, 

and  consequently,  q0  =  0. 

In  order  to  have  r0  -  0,  we  must  set  ~  0*. 
Then,  for  a  perfectly  regular  departure: 
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{  , 


6C  3  Vo  =  Yl  =  Po  =  9o  =  ro  “  0. 


V.  7 


Now,  let  us  find  out  whether  M/6  and  co,/6  have  limits  for 
t  =  0.  b 

If  M/6  has  a  limit,  a>,/6  will  have  a  limit,  and  it  will 
be  the  same  limit.  The  expression  for  M/6  is: 

M  =  _  q$  slnJL  +  3^  cos  0  cos  y  -  sin  9. 

The  last  term  vanishes  at  the  origin,  as  %  does. 

The  first  term  vanishes  at  the  origin,  as  yj,  does.  To  see 
this,  it  is  sufficient  to  apply  L* Hospitals  rule  to  the  ratio 

.s.j-.njL  an(j  utilize  formula  (15) . 

Now,  let  us  study  the  second  term.  By  (13) » 

-u  =  R'  (k  -  1)  I3-  , 

while 

6X  =  -  6  sin  y  . 

As  sin  y  vanishes  at  the  origin,  so  does  the  second  term. 

So,  in  the  case  of  perfectly  regular  departure. 


lim  M  =  lim  _  q 
t=o  ^  t=o  ~5 


(16) 


12  z 


Now,  we  must  calculate  the  derivative  -j^r  -  p? . 

Equation  (9)  gives: 

pt  =  5«  +  ©ti  cos  ^  +  Tjtt  cos  9  sin  y  -  7]f  9»  sin  9  sin  y+  oo^  y' 
Let: 

$  f  9"  cos  y+  >1"  cos  9  sin  y  -  j/Q*  sin  9  sin  y  ,  (17) 

6  =  Q«  cos  y  +  Ji"  cos  9  sin  y  .  (18) 

The  derivative  of  p  will  then  have  the  form: 

p»  s  511  +  $  +  *  yr  .  (19) 

4 
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At  the  origin: 

lim  |4  =  rim  y » 
t»o  t=o 

So,  the  terms  ^4  and 
factors  which  vanish  at  the  origin.  Therefore,  if  t  is  extreme¬ 
ly  small, 

„  S"  +  8  (  wlth  t  ■»  0>  (20) 


=  lim  =  lim  sin  y  =  0. 
t=o  ®  t=o 

T  9f  sin  9  sin  y  each  contain  two 
6 


and  the  closer  t  is  to  zero,  the  more  satisfactory  the  approxi¬ 
mation  will  be. 

6.  USE  OF  THE  MOMENT  OF  MOMENTUM  THEOREM.  Let  OX,  OY, 

('  ^  OZ  be  absolute  axes,  and  GX',  GY*,  GZ’  be  axes  respectively 

parallel  to  the  former.  Here,  the  absolute  axes  are  related 
to  Earthj  this  is  a  sufficient  approximation.  The  motion  with 
respect  to  the  trihedron  of  reference  GX'Y'Z'  is,  by  definition, 
the  '’motion  with  respect  to  the  center  of  gravity".  With  this 
’  new  trihedron  of  reference,  we  can  use  the  fundamental  moment 

of  momentum  theorem,  without  modifying  the  external  forces:  The 
derivative  of  the  vectorial  moment  of  momentum  (with  respect  to 
t)  is  equipollent  to  the  vectorial  sum  of  the  moments  (with 
respect  to  G)  of  the  external  forces*. 

We  do  not  have  to  bother  ourselves  with  fixing  the  absolute 
trihedron  OXYZ;  we  are  going  to  take  the  projections  on  the  axes 
of  the  free  trihedron,  which  has  been  defined. 

The  components  of  the  moment  of  momentum  on  Gx,  Gy,  Gz 

are: 

BP.*  AP 


The  components  of  the  relative  rate  of  change  (with  respect 
to  Gxyz)  are: 

*PT,  flq1,  Ap* • 

/  r  / 

•  Sea  e.g. ,  Appell's  Tralte  de  Uacanlque.  or  else  Painleve'a  Coura  de_ 
Mecanlque.  Gauthier-Villara,  1930. 
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The  components  of  the  translational  rate  of  change  are  the 
minors  of  the  matrix: 


,P  <3  r 

I  flp  Bq  Ap  I  * 

i  *  e  •  , 

dqp  -  Sqr,  -( A  pp  -  Apr) ,  flpq  -  5pq. 

We  have,  then,  the  three  equations: 

B  ’+  iiqp  -  flqr  = 

^dt“i4PP+^Pr=  (21) 


I  & 

T  rit. 


M.  Esclangon,  who  utilized  the  moment  of  momentum  theorem 
in  another  form,  fortunately  discovered  that  an  algebraic  equa¬ 
tion  of  the  second  degree  in  y1  can  be  rigorously  deduced  from 
these  equations.  The  condition  of  reality  will  give  a  lower 
limit  of  2. 


Before  seeking  this  equation  of  the  second  degree,  I  shall 
make  a  remark.  The  spin  of  the  projectile  has  a  negligible  varia 
tion  along  the  trajectory.  So,  the  role  of  M  is  almost  negli¬ 
gible.  To  simplify  the  explanation,  I  assume  that 

r 

M  »  0,  whence  p  =  const. 

*  ■  .  .  z  ,  .  r  ■  ,  .  . 

In  the  case  of  perfectly  regular  departure, 

Po  *  <Pl  =  2* 

In  the  contrary  case, 

p„  =  2  Q  . 


(22) 

(23) 


In  fact,  the  Initial  values  60,  y’,  7)f  are  extremely  small 
if  the  projectile  is  well  designed  and  if  the  tube  is  not 
eroded. 
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Note  that,  if  is  not  completely  negligible,  it  is  suffi- 

cient  to  add  to  the  expression  for  p  a  function  D(t),  which 
vanishes  at  the  origin,  decreases,  and  is  extremely  small  if 
the  projectile  is  normal. 

The  first  of  equations  (21) : 

5p*  +A  qp  —  #qr  —  *x>  (24) 

would  be  an  absolutely  rigorous  equation  if  the  parameters  k 

and  l  ,  which  figure  in  M  .  were  known.  These  parameters 

^ . 

certainly  vary  along  the  trajectory;  but  they  vary  little  if 
the  stability  is  suitable,  so  that  the  drift  is  regular  and 
feeble.  If  we  assume  that  k  increases  between  the  origin  and 
the  dangerous  zone,  and  that  it  can  take  values  between  1  and 
3  or  4 ,  approximately,  this  is  an  hypothesis  which  appears  ac¬ 
ceptable.  Furthermore,  it  is  possible  that  the  variations  of 
l  are  in  the  sense  opposite  to  those  of  k. 

In  the  case  of  flat  fire,  it  is  extremely  probable  that 
the  variations  of  k l  are  quite  small.  It  is  difficult  to 
give  general  indications  from  this  point  of  view  in  the  case 
of  very  high  angle  fire. 

We  could  study  equation  (24)  with  great  precision  by  keep¬ 
ing  sin  6  and  cos  6,  by  not  assuming  any  approximation.  How¬ 
ever,  that  would  be  difficult.  So,  I  shall  make  certain  approx¬ 
imations  which  will  not  alter  the  character  of  the  equation  at 
all. 

Let  us  consider,  e.g.,  the  derivatives  9T  and  ft . 

At  the  origin,  9for  V:  will  be  very  small  if  V0  is  very 
large. 

Later,  9’  or  t  f  will  remain  very  small  if  V  remains 
quite  large.  I  shall  always  suppose  that  V  remains  quite 
large,  so  I  exclude  not  only  the  case:  a=  90°,  but  also 
the  case:  a  *90°.  I  do  not  study  almost  vertical  fire  because, 
near  the  summit,  the  minimum  velocity  would  be  quite  close  to 
zero.  The  phenomena  then  become  entirely  different  from  those 
which  we  are  examining. 

Now,  let  us  consider  the  derivative  rj1  . 

If  6°  -  0,  or  6°  =  0,  we  shall  have:  pj  =0. 

i 

In  the  case  of  perfectly  regular  departure,  r)^  =  0. 

In  the  contrary  case,  will  be  very  small  at  the  origin, 
as  is.  Later,  will  remain  small  if  Vm  remains  quite  large 
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The  variations  of  k  and  l  play  a  role  here.  Let  us  repeat 
that  this  role  is  a  minimum  in  the  case  of  flat  fire  if  the 
angle  of  projection  is  large. 

As  long  as  V  remains  large,  our  ignorance  on  the  sub¬ 
ject  of  k  and  l  m  probably  gives  us  no  great  embarassment . 

When  I  say  that  the  minimum  of  V  must  be  quite  large, 
this  statement  is  not  precise  enough.  Suppose  V0  very  large 
and  a  quite  near  90°.  The  summit  S  will  have  a  high  alti¬ 
tude  and,  consequently,  at  the  point  of  minimum  velocity, 
the  ballistic  coefficient  c  will  be  very  much  diminished  with 
respect  to  its  initial  value  c0.  There  results  a  great  dimi¬ 
nution  in  the  value  of  the  stability  function  Q.  So,  the  nu¬ 
merical  value  of  the  mean  yaw  might  become  too  large,  with  regard 
to  the  general  conditions  of  stability.  I  shall  say,  then: 

V  must  always  remain  quite  large,  and  all  the  larger,  the 
hxgher  the  corresponding  point  of  the  trajectory  is. 

So,  I  always  suppose  that  the  initial  velocity  V0  is  very 
large,  that  the  minimum  velocity  is  not  small,  and  that  the 
angle  of  projection  a  is  different  from  90°.  Otherwise,  the 
stability  is  doubtful  a  priori;  absolutely  nothing  would  be 
known  about  the  parameters  k  andr  l ;  certain  hypotheses  such  as: 

6  »  t  ,  ©»  w  ,  9n  «  Tn, 

might  have  no  foundation.  Besides,  there  would  be  no  sense 
in  writing: 

sin  6  *  6,  and  cos  6  *1. 

Without  actually  insisting  on  these  questions  of  approxi¬ 
mation,  I  must  underline  their  importance. 

Let  us  analyze  all  the  terms  of  equation  (2 4) • 

Neglecting  the  product  we  have: 


pi  «,  fin  +  @  +  yt . 

Neglecting  products  such  as;  9'68,  qr6z,  9?26,  q,26, 

11*9' 6,  etc.,  we  have: 

qr  *  6  YtS  +  (2co^6  +  co^) . 

Therefore, 

^  ;  p»  -  qr  *  6”  +  0  ,  -  6  f  ’ 3  -  2C026H1  •  (25) 
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On  the  other  hand,  if  we  suppose  M  =  0, 

z 


P  =  Po  =  S2  x  - 


The  very  approximate  expression  for  q  is: 
q  -w  6  y»  +  M. 


So: 

qp  *  0^(6  y’  +  co^6  +  <o^)  .  (26) 

The  first  member  of  equation  (24)  will  be,  with  a  good 
approximation  if  we  remain  within  the  limits  previously  in¬ 
dicated: 


fl&j.(6  y*  +  +  ®^)  +B:  (6”  +©  -  8  f 2  -  2oo^6  y.1 )  • 

However,  we  know,  a.  priori ,  that  2  !  must  be  a  large  num¬ 
ber  and,  on  the  other  hand ,  B /A  is  certainly  less  than  8  or 
10,  for  example. 

Besides,  [co^j  remains  very  small. 

So,  we  can  certainly  neglect  the  last  term,  for  it  is 
always  extremely  small  with  respect  to  the  first.  So,  for 
the  primitive  equation  (24) »  we  can  substitute  the  following 
approximate  equation; 

'  -i  • 

-  :56y»8  +  A  0^  yi  +  B( 6"  +0  )  +  A  -  Rk  l  6  =  0.  (27) 


Roughly,  the  approximation  is  good  if  6  remains  small, 
and  the  better,  the  closer  6  is  to  zero. 

The  condition  of  reality  of  y1  will  be; 

(A  Gi6)2  +  -4  5  6  [5  (8«  +  0  )  +AQ  xM  -  Rk  l  6]  =  0.  (28) 

I  recall  that: 
j  R  =  mSr  =  mcF(v) , 
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o 


©  =  9”  cos  \j/  +  »]"  cos  9  siny  , 

and 

o>^  =  -  9*  sin  y.  +  *)*  cos  9  cos  y  . 

This  is  N.  Ssclangont s  criterion,  inutatis  mutandis  (loc. 
cit.,  p.  780). 

I  shall  put  it  in  the  following  form: 

SS2  +  u  Bz  6  (6»  +  e  )  +  J^ABQ^n  =  0,  (29) 


setting 

S=(lQ1)2-45Rki. 


Since  the  yaw  5  is  positive,  the  criterion  (29)  can  be 
put  in  the  form: 


SS  +  4  b  2(6«  +6  )  +  4a  BQ.  iM  =  0. 

(30) 

M.  de  Sparre,  in  his  aforementioned  Report  of 
and  295) ,  gave  two  successive  criteria: 

1904  (pp.  293 

( A  £2  ) 2  >  4  B  (Rk  i  )  0 ,  .  f  ' 

(31) 

or 

S0  >  0,  with  £2j_  =  G  , 

and 

{A  £2  )2  >'  y-  B  (Rk  i  )0\ 

(32) 

The  first  criterion  is  deduced  from  the  equations  of 
gyroscopic  motion,  derived  in  the  Report  of  1904,  and  the 
second  criterion  is  deduced  from  very  ingenious  considerations 
of  the  physical  character  of  the  air  resistance. 

I  shall  study  criterion  (30)  at  the  origin,  l.e.,  for 
t  «  0,  As  a  result  of  this,  if  the  departure  is  perfectly 
regular,  the  second  criterion  is  applicable  again,  for  t  *  0. 


( 
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As  a  result  of  this.  If  the  departure  Is  perfectly  regular,  the 
second  criterion  Is  applicable  again,  for  t  v  0.  However,  if  the 
departure  has  the  least  irregularity,  it  must  be  foreseen  that 
the  second  criterion  will  be  insufficient  for  t  *  0.  M.  de  Sparre 
always  implicitly  supposed  that  the  departure  was  perfectly  regular 

After  the  criterion  has  been  studied  at  the  origin,  it  must  be 
studied  for  any  values  of  t  whatever;  but  actually,  I  put  this  ques 
tion  aside. 

M.  Esclangon  made  a  penetrating  study  of  this  question  by  dis¬ 
tinguishing  _2onas__of_stability  and  zones  of  instability  on  a  tra¬ 
jectory. 

7.  INTEGRATION  OF  A  DIFFERENTIAL  SYSTEM.  I  suppose  that 
the  departure  is  perfectly  regular  and  that  M  =  0.  Under  these 
conditions,  equations  (2l)  become: 


B  p’  +  q  04  £2  -  Jr)  =  #x, 
B  q'  -  P  (A  a  -  B  r)  =  Hy, 


(33) 


with 


P  -  Po  =  a  • 

Besides, 

r0  -  0. 

So,  the  term  #  r  is  negligible  with  respect  to  A  2  at  the  origin, 
i.e.,  when  t  *  0,  and  the  above  equations  can  be  replaced  by  the 
following  with  an  excellent  approximation: 


B  p»  +  q4£2=^x, 

jq’  -  VAQ,  =  My»  (34) 

These  are  the  differential  equations  which  M.  Sugot  very 
fortunately  used  to  derive  the  approximate  equations  of  gyroscopic 
motion.  His  equations  are  similar  to  those  of  M.  Esclangon  and 
the  approximations  are  the  same.  I  am  going  to  study  equations 
(34) t  so  as  to  be  able  to  discuss  criterion  (30) . 

Let  a)  =  aQ*/b  ;  (0  will  be  a  large  number,  like  Q  *  I  suppose, 
in  fact,  that  the  ratio  B  / A  is  less  than  7  or  8',  for  example. 


We  shall  write: 


.  V- 


+  ooq  =  <oP6  -  fXj 
^  -  wp  =  <dA  8  -  fa. 


(35) 


Although  it  appears  a  bit  paradoxical  at  first  glance,  I 
am  going  to  integrate  equations  (35)  as  if  fx  and  f2  were  knomi 
functions  of  t.  These  functions  are  continuous  and  differentiable, 
since  we  suppose  essentially  that  the  derivatives  5*  and  6”  exist, 
in  all  that  precedes. 


Note  that  the  ratio  fg/fj.  -  A  /P  is  supposed  constant  and 
very  small  for  t  *  0.  We  assumed  that  this  ratio  can  have  a  value 
such  as  2 %  or  3%  at  the  origin. 

The  classic  method  of  ''variation  of  constants”  immediately 
gives  the  solution: 


p  =  X  cos  arc  +  p.  sin  cot, 
q  =  Xsin  cot  -  p  cos  cot. 


(36) 


o  Whence: 

Xf  =  fi  cos  cot  +  f2  sin  cot, 
(a>  =  fi  sin  cot  -  f2  cos  cot. 


(37) 


To  have  p0  =  q0  s  fe  must  take  Xo=|ao  =  0.  It  will  be 
convenient  to  introduce  the  following  symbols:  , 


f i /co  —  F  x ,  f  2  /go  —  Fg. 


It  should  be  remarked  that,  always, 
pa  +  q2  =  X  2  +  fi  2 . 

I  shall  have  to  be  occupied  with  the  case  where  Fx  and  F2 
are  linear  functions  of  t: 

Fx  (t)  =  Hxt,  Fz(t)  =  Hat,  H2/Hj.  =  A/P.  ,  (38) 

It  will ‘be  convenient  to  introduce  the  auxiliary  variable: 
T  55  •  cot . 
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Immediately: 

T  T 

co  X  =  Hx  /  T  cos  T  dT  +  Hg  /  T  sin  T  dT, 
o  o 


co  X  =  HX(T  sin  T  +  cos  T  -  1)  +  Ha  (-T  cos  T  +  sin  T),  (39) 

and,  similarly,  • 


co  jk  =  Hi(-  T  cos  T  +  sin  T)  -  Hg(T  sin  T  +  cos  T  -  1) .  (40) 

The  curves  of  X  and  p  can  be  constructed  according  to  the 
formulas,  remarking  that  the  extremes  depend  uniquely  on  the 
ratio  A/P,  from  the  point  of  view  of  their  position. 

On  the  ot*  3r  hand,  the  terms  in  H2  will  generally  be  negli¬ 
gible  with  respect  to  the  terms  In  Hx,  because  the  constant  Ha 
will  be  a  very  small  number,  and  because  the  ratio  A /P  is  small 
with  respect  to  unity. 

I  adopt  the  following  notation: 


for 


or 


x  =  x0,  Xj,  x2,  x, ,  x^,, 


T  =  0,  |  !*,  2*, 


_  n  *  &  3*  2* 

x  ~  u»  2co>  ci)'  2co*  co  ‘ 


Hence: 


a)  \1  =  Hx  (f  -  1)  +  H2,  cop  =  H  -  (§  -  1)H2, 


coX  2  —  -  2Hi  +  #H2,  cop  2  ~  +  2H2 , 


«  =  -  (3§  +  1)H,  -  H2,  ton,  =  -  Hx  +  0%  +  l)Ha, 


co  X  ^  —  —  2*H2 , 


co  p  .  =  -  2^Hi . 
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A 


Consequently: 

Pi  =  Pi“  fs 

p3  =  "  p3“  »L’ 


Pa  -  -  Xe 

p4  “  x4  “  z*t’ 


(41) 


Qi 


C3 


q»  =  Ma  * 


(3|  +  «5*. 


K,  =  2#— 1 

r4  .CO 


(42) 


The  curves  of  p  and  q  can  be  constructed,  remarking  that 
the  extremes  of  p  are  given  by  the  equation: 

q(t)  =  Fx(t),  (43) 

and  the  extremes  of  q  by  the  equation: 

")  P(t)  =  -  F2(t),  (44) 

so  that  the  extremes  of  p  correspond  to  the  following  values: 

T  =  0,  T  *  it,  and  T  =  2*. 

The  extremes  of  q  correspond  to  the  values: 

T  =  0,  T  =  2#. 

Figures  6  and  7  represent  the  curves  of  X  f  M ,  p,  and  q 
in  the  interval  from  T  =  0  to  T  =  2*. 

Besides,  we  must  note  this:  whatever  the  functions  P(t) 
and  6(t)  are,  if  they  are  finite  and  continuous  in  the  neighbor 
hood  of  the  point  t  =  0, 

lira  £  =  lim  £  =  0. 
t=o  6  t=o  6 

In  fact,  let  us  return  to  equations  (36)  and  (37),  and 
suppose: 

X  O  Po  —  Po  ~  Qo  ~  0  • 


t 
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The  functions  and  fz  contain  6  as  a  factor.  It  will  be 
sufficient  to  apply  the  theorem  of  the  mean  in  the  interval  from 
t  =  0  to  t  =  e,  and  to  make  e  approach  zero  as  a  limit. 

Remark.  I  must  make  a  remark.  In  this  theory,  there  exists 
an  appearance  of  a  vicious  circle,  for  I  seek  the  value  of  2  and 
at  the  same  time,  I  sometimes  rest  on-  the  fact  that  £2  will  be  a 
large  number.  However,  this  is  only  an  appearance  of  a  vicious 
circle.  In  fact,  2  is  not  a  parameter  that  is  susceptible  to 
vary  from  zero  to  infinity.  With  a  given  caliber  and  initial  vel 
ocity,  we  are  certain  that  2  will  vary  between  two  rather  close 
limits  2»  and  2*.  That  is  why  we  can  regard  2,  a  priori ,  as  a 
large  number  with  respect  to  other  elements  which  come  into  play. 

Finally,  when  2  will  have  been  determined  according  to  the 
criterion,  we  shall  see  well  whether  the  method  followed  is  good, 
mediocre,  or  bad.  In  the  last  case,  the  method  will  have  to  be 
perfected. 

For  example,  the  major  calibers  might  require  particular 
studies;  but  I  am  not  actually  making  a  complete  discussion  of 
all the  cases. 

8.  STUDY  OF  THE  YAW  AT  THE  ORIGIN.  I  suppose  that  the 
departure  Is  perfectly  regular;  I  am  going  to  make  the  approx¬ 
imate  calculation  of  5  at  the  origin,  i.e.,  when  t  varies  from 
0  to  2«/co. 

I  shall  make  a  first  approximation  by  starting  from  formulas 
(9)  and  (15) ,  and  setting  - 

6  *1  t»  j  t.  (45) 


That  will  permit  the  knowledge  of  p,  with  a  great  approxi¬ 
mation,  according  to  the  preceding  study  of  p  and  q. 

*  »  . 

Knowing  p,  I  shall  make  a  second  approximation  of  6  by 
setting 


6»  *|  t  »|  +  p. 

C 


(46) 


That  will  give  the  form  of  the  curve  S(t)  at  the  origin. 

Note  that  t  varies  here  in  an  extremely  small  interval  if 
(o  is  a  large  number.  As  a  result  of  this,  we  know,  we  can  regard 
P(t)  and  T»(t)  as  constants  and  equal  respectively  to  P<,  and  t» 
in  the  Interval  from  T  =  0  to  T  =  2*.  ° 


So,  the  first  approximation  will  give 


P5«Pe|t;i  t.  (47) 


If  we  refer  to  the  previous  notations,  we  must  set,  here: 
Hi=P„|T>|  ,  or  I1  =  ;i.  (48) 


I  suppose  that  the  ratio  P0/a>  is  very  small.  (Besides, 
t  2>  is  very  small  if  V0  is  very  large)  .  Consequently,  the  ordi- 
hates  of  the  curve  of  p  (Fig.  7)  are  always  very  small  with  res¬ 
pect  toT £ .  The  curve  which  represents  6’(t)  will  he,  with  a 
very  great  approximation,  the  curve  of  p,  hut  referred  to  an  axis 
0rX’,  parallel  to  OX,  the  distance  00*  being  equal  to  |t  »  i(Fig.  8). 

The  function  5’  will  have  a  maximum  for  T  =  tc  +  h,  where 
h  w  0,  and  a  minimum  for  T  =  2*. 

So,  the  curve  of  5  will  have  a  point  of  inflection  for  T  = 

#  +  h,  and  another  for  T  =2#. 

Here,  then,  is  an  essential  fact:  in  the  interval  from 
T  =  0toT  =  #  +  h,  5n>  0,  and,  in  the  interval  from  it  +  h  to 
2*,  5"  <  0. 

When  t  varies  from  0  to  about  tt/oi,  the  curve  5(t)  turns  its 
concavity  upwards,  and  when  t  varies  from  about  «/co  to  2tt/ co,  the 
curve  6(t)  turns  its  concavity  downwards;  this  is  an  Important 
fact. 


9.  STUDY  OF  THE  CRITERION  AT  THE  ORIGIN.  DEPARTURE  PER¬ 
FECTLY  REGULAR.  I  suppose  that  the  departure  is  perfectly  regu¬ 
lar.  So,  it  will  be  necessary  to  substitute  Qfor  ^  in  the 
inequalities  (29)  and  (30) . 

Besides,  I  actually  study  the  criterion  only  at  the  origin, 
i.e»,  for  t  *  0. 

By  the  relation  (16) , 

lim  M  =  0, 
t=o  5 

and  M/6  is  negligible  for  t  0  in  the  inequality  (30). 

If  there  is,  for  t  *  0,  a  double  inequality  of  the  type: 


0< 


(49) 


<  A, 


where  A  is  positive,  then,  at  the  origin,  we  shall  be  able  to 
adopt  the  criterion  (31).  Moreover,  it  must  be  remarked  that 
this  criterion  is  obtained  immediately  if  it  is  assumed  that, 
at  the  origin,  the  projectile  is  similar  to  the  top  —  a  rather 
fragile  intuition,  when  it  concerns  so  delicate  a  discussion*. 


However,  I  have  Just  shown  that  5”  changes  sign  in  the 
immediate  neighborhood  of  the  origin.  As  a  result  of  this,  the 
possibility  of  this  relation  (49)  is  not  probable,  and  I  am  go¬ 
ing  to  study  carefully  the  values  of  6”  +  fo  for  £  «  q. 

'6 

With  reference  to  formula  (20): 


£'  *  5"  +  e 
6  6  ' 


for  t  *  0} 


it  will  be  sufficient  to  study  the  values  of  p’/S  for  t *  0. 
-With  a  sufficient  approximation, 

pt  =  oj(po6  -  q),  (50) 


f  =  a)(P„  -  J)  .  (51) 

So,  it  is  necessary  to  study  the  variations  of  q/6  for 
t  ^  0 » 

In  the  interval  from  T=0toT=«+h,  p*>  0.  However, 
in  the  interval  from  co  +  h  to  2#,  pf<  0. 

So,  it  is  necessary  to  know  the  maximum  of  the  modulus  of 
.  (Pe6  -  q)  in  this  interval;  it  is  this  maximum  which  intervenes 
in  the  criterion. 

The  following  values  are  very  approximate: 

for  T  =  *,  6  =  *  I  ,  q  =  £  |  t’|  p„; 
for  T  =  2»,  8  =  ,  q=f*l*;|P.. 


*  On  the  subject  of  the  elementary  theory  of  tho  top,  9ee  e.g. ,  Beghin, 
Statics  and  Dynamics,  Vol.  II,  collection  Armand  Colin,  p.  168. 
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So,  at  these  two  points,  the  value  of  q/5  is  the  same,  with 
a  great  approximation. 

The  following  values  are  also  very  approximate: 

for  T  =  §*,  6  =  |*  I  t’|  ,  q  =  +  1)  I  T*  |  P„. 

So,,  at  this  point,  the  value  of  q/6  is  greater  than  at  the 
points  *  and  2if.  If  the  number  P0/co  is  very  small  (which  can  be 
verified  when  the  value  of  Q  is  fixed) ,  all  these  approximations 
are  satisfactory:  but  they  are  approximations;  so  it  is  useless 
to  seek  the  exact  maximum  of  q/6.  We  have  immediately  a  superior 
limit  of  this  maximum,  by  taking  the  largest  value  of  q  in  the 
interval,  and  the  smallest  value  of  6.  That  gives: 


sup  lira  q/6  =  2P0,  (52) 

sup  lim  (q/6  -  P0)  =  P0  (53) 

in  the  interval  from  T  3  k  +  h  to  T  =  2#,  or,  finally, 

sup  lim  |  pT/M  -  oiP0»  (54) 

sup  lira  I  1~  °>po  (55) 


The  inequality  (30)  can  then  be  put  in  the  following  form 
for  t  *  0: 

S  =  4  fl2  (oP„,  (56) 


where 

S-(4fi)z-4  B  (Rk  l  )  0.  ; 


i.e., 

(A  Q ) 3  =8S  (RkZ  )0, 


(57) 


for 

P  =  Rk* 

A  Q  *  . 

The  coefficient  8  in  (57),  which  I  obtain  thus,  is  rather 
high.  M.  de  Sparre’s  coefficient  in  (32)  is  16/3.  However,  I 


immediately  find  M.  de  Sparre’s  coefficient  again  intuitively. 
According  to  the  preceding,  the  maximum  of  q/6  will  occur 
approximately  at  T  =  3«/2.  At  this  point,  with  a  great  approx¬ 
imation, 

q3  =  <■¥  +  15  5|T°I  >  63  =  25  1  ^  I  > 


3*/2  +  1 
3k/2 


y 


If  I  substitute  1/3  for  2/3#  >  I  obtain  M.  de  Sparre’s 
coefficient.  As  this  theory  has  no  point  of  contact  with 
M.  de  Sparre’s,  I  consider  the  concordance  of  figures  as  be¬ 
ing  very  satisfactory.  However,  this  accord  exists  only  in 
the  exceptional  case  of  perfectly  regular  departure. 


10.  INTEGRATION  OF  A  DIFFERENTIAL  SYSTEM  NEAR  THE  PRECED¬ 
ING  SYSTEM.  We  shall  now  have  to  consider  a  non- van! s hlng 
initial  yaw.  As  a  result  of  this,  it  will  be  opportune  to 
integrate  the  system  (35)  by  taking  f-L  and  fz  as  constants. 

I  obtain  the  expression  of  these  constants  by  taking  the 
values  of  P  and  5  at  the  instant  t  =  «/co,  for  my  study  will 
be  limited  to  the  Interval  of  time  from  t  =  0  to  t  =  2#/a>. 

By  thus  taking  constant  approximate  values,  I  introduce  an 
artificial  periodicity;  but  that  offers  no  inconvenience  at 
all,  since  this  study  is  actually  limited  to  the  interval  of 
time  indicated  above. 


Let 


Hence, 


fx  =  coAj  and  =  00A3. 


A2/A1  =  A/P, 

and  we  know  that  this  ratio  is  very  small.  For  example,  at 
the  origin,  its  value  might  be  1 %  or  2%  or 

First,  suppose: 

-Ho  —  Po  —  Qo  —  0. 


Integration  of  equations  (37)  give's  immediately: 


(58) 


X  =  Ax  sin  cot  +  A2  (1  -  cos  cot) , 
p  =  AiCl  -  cos  cot)  -  A2  sin  cot. 


It  is  convenient  to  introduce  the  variable  T  =  cot. 

It  is  seen  immediately  that  X  has  a  maximun  near  Ai  for 
T  *  */2,  and  a  minimum  near-  -Ax  for  T  *  3 «/2,  since  A2  is 
negligible  with  respect  to  Ax. 

Similarly,  p  has  a  maximum  near  2Aa  for  T  w  rt,  and  a  very 
small  minimum  for  T  *  0.  Also,  it  is  evident  that  the  curves 
of  X  and  p  are  very  close  to  each  other,  and  that  the  curves  of 
ja  and  q  are  very  close  to  each  other.  That  originates,  par¬ 
ticularly,  in  the  following  fact:  - 

p8  +  q8  =  X2  +  p8  =  4(A2  +  A2)  sin  2  (cot/2) ,  (59) 

'*■  1  2 

so  that  the  points  (X,  p )  on  the  one  hand,  and  (p,q)  on  the  other, 
describe  two  circles  of  the  same  diameter,  passing  thru  the  ori¬ 
gin  (Fig.  10) .  If  we  consider  the  tangents  at  the  origin,  these 
two  tangents  are  symmetric  with  respect  to  Ox,  and  very  slightly 
inclined  to  Ox.  So,  the  centers  Bx  and  B2  of  the  two  circles 
are  symmetric  with  respect  to  Oy,  and  are  very  close  to  it;  they 
would  be  confused  if  the  ratio  Az/Ax  were  zero.  The  symmetry  of 
the  tangents  with  respect  to  Ox  results  immediately  from  the  follow¬ 
ing  expressions:  .. 

pf  -  coAi,  qf  =  C0A2,  X 1  =  coAx,  p  ’  =  -  0A2,' 


*’/*•  -  -  q’/p* 

#  O  O  a 


It  must  be  noted  that  the  forms  of  the  four  curves  (Fig.  9) 
are  determined  by  the  value  of  the  ratio  Az/Alf  and  this  is  what 
makes  this  approximation  interesting. 

Now,  suppose: 

Po  /  0,  q  /  0. 


Since 


we  will  have: 


Po  -  Xc  and  q0  -  -  p0j 
*  o  /  0  and  0  /  0 . 
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Let  us  keep  the  notation  X ,  p  ,  p,  q  for  the  expressions 
which  have  been  calculated  and  which  vanish  at  the  origin.  We 
must  now  replace  X  and  p  by  X+  XG  and  p  +  p0. 

The  new  expression  of  p  will  be  As 

A  =  p  +  X 0  cos  cot  +  p  0  sin  cot. 

The  new  expression  of  q  will  be  B: 

B  -  q  +■  X0  sin  cot  -  p0  cos  cot. 

What  has  been  done  permits  us  to  integrate  the  system  (35) 
in  somewhat  more  general  conditions. 

Let  (px,  pi)  be  the  solution,  which  vanishes  at  the  origin, 
for  the  case: 

fz  =  coA  i  ,  fa  =  C0A2 . 

Let  (p2,  qz)  be  the  solution,  which  vanishes  at  the  origin, 
for  the  case: 

=  coHit,  fz  =  coH2t; 

(pi  +  Pa>  Qi  +  qz)  will  be  the  solution,  which  vanishes  at  the 
origin,  for  the  following  case: 

fi  “  co(Ai  +  Hit)  ,  fz  =  co(Az  +  Hzt)  . 

Furthermore,  to  have  a  solution  that  passes  thru  the  point 
(Po,  qD)  at  the  origin,  it  is  sufficient  to  add  the  periodic 
functions:  ' 

Wx  -  cos  cat  +  40  sin  cat,  (60) 


Ws  =  X0  sin  cat  -  p0  cos  cat,  (61) 

the  solution  of 
and  the  general 
the  origin,  by  A 


11.  COMMENTARY  on  IRREGULAR  DEPARTURES.  It  must  be  fore¬ 
seen  that  the  departure  could  not  be  perfectly  regular.  It  cer¬ 
tainly  depends  on  us  to  utilise  only  well  inspected  shells  and  un 
eroded  tubes,  but  other  facts  do  not  depend  on  us.  Violent  re¬ 
actions,  which  make  the  trunnions  of  the  gun  vibrate,  can  be 


To  avoid  all  confusion,  we  could  designate 
(35) ,  which  vanishes  at  the  origin,  by  p  and  q, 
solution,  which  assumes  the  values  p0  and  qe  at 
and  B . 
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produced  at  the  contact  of  the  earth  and  the  carriage.  Further¬ 
more,  at  the  muzzle,  the  expansion  of  the  gases  and  their  in¬ 
flammation  cap.  generate  certain  vortices  and  cavitations  which 
act  on  the  projectile  as  deviating  forces.  From  the  analytic 
point  of  view,  we  say  that  an  irregular  departure  corresponds  to 
a  non-vanishing  initial  value  of  the  yaw,  and  an  initial  value 
d&  t 

of  -pj:  different  from  the  regular  initial  value  |  td|  . 

However,  a  distinction  must  be  made. 

If  very  violent  perturbations  existed  at  the  departure, 
we  would  then  have  to  study  the  primitive  equations  (21),  and 
not  the  equations  such  as  (27) ,  deduced  from  the  primitive 
equations  by  simplifications  which  are  permitted  by  reason  of 
a  certain  regularity  of  the  physical  phenomenon,  of  which  we 
have  a  well  founded  intuition.  The  study  of  very  violent  per¬ 
turbations  is  not  a  question  which  properly  concerns  ballistics: 

I  do  not  speak  of  it* 

I  am  going,  to  consider  exclusively  the  following  conditions, 
which  seem  to  me  normal  if  the  initial  perturbation  is  small. 

I  assume  the  existence  of  an  extremely  small  initial  value 
60,  say  of  the  order  of  a  thousandth  of  a  degree,  this  figure 
being  only  an  indication. 

Furthermore,  I  assume  the  existence  of  a  modification  of 
the  regular  value  of  which  is  |  t»|,  a  small  modification 

compared  to  the  regular  value  [ |  ♦  ° 

O 

I  shall  remain  almost  in  the  preceding  frame  of  study, 
and  shall  say  that  the  departure  is  almost  perfect. 

O 

To  simplify  the  facts,  I  suppose  y0  =  0,  i. e. ,  6  -  0. 
Whence,  60  =  &£,  while  60  =  0  in  the  case  of  perfectly  regular 
departure. 

As  a  result  of  this,  pi  =  0. 

Consequently,  as  in  the  case  of  perfectly  regular  departure: 


0.0  o  o  o 


M°  =  N°  =  0. 


Let  us  consider  the  expressions  for  p,  q,  p.  We  have: 

i  * 

Po  =  &!>  +  9’  “  &'  +  •  (62) 

o  o  o 


I  assume  again  that  0’  =  t1,  because  60  is  extremely  small. 


Then: 


vO 

0 

II 

Vi  sin  6o  WVi6o» 

•  (63) 

T0  = 

yi  cos  50  ft*  yk, 

« 

(64) 

Po  = 

Q  +  T0. 

(65) 

Let  us  return  to  the  fundamental  system  (21),  supposing 
M„  =  0,  for  greater  simplicity. 

Zi 

The  third  equation  will  give: 

(66) 

(67) 


p  -  Po  r  ®  ^ 
P  *  Q  +  vl. 


I  shall  suppose  essentially  that  the  ratio  |r0|/co  is  small; 
consequently,  the  ratio  jr0(/£2  will  be  a  very  small  number.  We 
shall  then  have  a  sufficient  approximation  if  we  substitute  Q 
for  p  in  equations  (21),  and  4p  for  4 p  -  B r,  for  t  *0.  So,  vje 
again  find  equations  (34),  but  with  not  quite  so  good  an  approxi¬ 
mation  as  in  the  case  of  perfectly  regular  departure;  for,  in 
the  latter  case,  rc  =  0. 


I  have  thus  substituted  the  following  for  the  system  (21): 


ft  +  tD<1  =  V 


(68) 


P  O  • 


The  system  (21)  is  not  an  ordinary  system  defining  p,  q,  p, 
because  6,  which  figures  in  the  second  members,  is  also  an  un¬ 
known  function  of  t.  However,  I  am  going  to  place  myself  in  par¬ 
ticular  conditions.  I  consider  only  the  very  small  interval  of 
time  from  t  =  0  to  t  =  2tf/a),  corresponding  to  a  very  small  zone 
Z0  on  the  trajectory.  In  this  zone  Zc,  I  shall  be  able  to  obtain 
a  very  satisfactory  approximation  of  6,  as  I  have  already  done, 
and  that  will  permit  me  to  regard  the  system  (21)  as  a  differential 
system  defining  p,  q,  and  p.  Of  course,  I  suppose  the  parameters 
k  and  J.  to  be  well  known  at  the  origin. 


furthermore,  considering  (21)  or  (68)  from  this  point  of 
view,  I  shall  be  able  to  choose  p0,  q0,  pQ  arbitrarily,  taking 
account  of  the  condition  which  I  have  imposed  on  myself  concern 
ing  the  limitation  of  |rJ/aj. 
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To  give  p0,  qG,  p0,  or  p0,  qoJ  re,  is  to  give  the  initial 
perturbation  of  the  free  trihedron. 


n 


I  essentially  suppose  that  this  perturbation  is  a  minimum. 
For  example,  the  ratio  I  p0/  ^£1  will  be  less  than  1/20  or  1/30. 

.  '  The  value  of  60,  e.g.,  could  be  less  than  1/10,000  or 

1/100,000,  which  corresponds  to  about  jq’qqq  or  ioo^OOO* 

I  shall  set  «{<{,  =  d  P0,  where  \a  I  can  be  comprised  between 
zero  and  a  number  quite  close  to  1,  for  example. 

If 

in  i=i/ioo, 

e.g.,  we  will  have 

I  Pol  <  1/2000. 

Then  take: 

60  =  1/60,000,  about  1°/100C)}  and  yj,  =  5; 
this  will  give: 

.  q0  =  1/12,000,  and  r0  *5. 


In  general,  the  figures  which  I  give  are  only  slightly  probable 
numbers.  The  last  ones,  in  particular,  constitute  only  a  vague  in¬ 
dication.  In  each  case,  it  will  be  necessary  to  examine  what  can 
be  a  normal  initial  perturbation.  Furthermore,  from  the  point  of 
view  of  approximations,  the  caliber  plays  an  important  role,  ac¬ 
cording  to  the  fundamental  formula: 


*  =  2Vj/a. 


The  design  of  the  shell  also  plays  its  role,  since 
o)  =  A  Q/B. 

We  have  previously  seen  the  role  which  the  ratio  P0/co  plays, 
and  I  supposed  that  this  ratio  is  small.  It  is 

P°  _  g  Pp 
o)  “  *nr, 

and,  if  a  designates  the  diameter. 


P0/Q  =  nacj 

so  that  P0/a)  is  generally  small,  because  the  ballistic  coefficient 
c  is  a  very  small  number.'  However,  it  will  be  necessary  to  examine 
all  the  possible  cases,  in  order  to  evaluate  the  coefficient  n. 

Since  the  departure  is  almost  perfect,  my  method  of  approxi¬ 
mation  will  be  modified  a  little.  Formula  (62)  gives 

=  I  n  1  +  Po. 


The  first  approximation  of  6f  will  be'  6*  =  5£,  which  gives 

6  =  60  +  (m  I  +  Po)  t.  (69) 

The  function  p,  or  A(t),  will  be  calculated,  taking  the  value 
pe  at  the  origin;  and  the  second  approximation  of  61  will  be 

6»  =  |  n  |  +  A(t).  (70) 

p,  or  A(t),  and  q,  or  B(t),  should  be  calculated  according  to 
the  system  (35),  in  which  the  second  members  comprise  a'  constant  and 
a  term  in  t,  according  to  the  approximate  expression  for  5,  given  by 
formula  (69) * 

We  shall  be  able  to  construct  the  curve  of  q,  or  B(t),  relative 
to  the  zone  Z0.  In  this  zone,  the  first  approximation  of  the  yaw  is 
represented  by  the  straight  line  0 ” e e 1 ,  whose  slope  is  6£;  and  0E  =  60. 

Draw  the  tangent  0>'E  to  the  curve  q;  the  ratio  of  the  slopes  of 
the  straight  lines  0»E  and  0ne gives  the  maximum  of  the  ratio  q/6 
(Fig.  11). 

The  constants  p0  and  q0  can  again  effect  a  raising  of  the  curve 
of  q  and,  consequently,  increase  the  maximum  of  q/6. 

It  will  not  be  useless  actually  to  consider  the  nature  of  the 
approximations  adopted  in  this  study.  Equation  (9)  is  an  exact  re¬ 
lation,  The  last  term  in  is  negligible  in  the  zone  Z0,  for  it 
contains  two  factors  which  vanish  at  the  origin.  We  can  substitute 
cos  y  0  for  cos  xf/„  ifv£  is  quite  small.  Besides,  here,  cos  V  o  =  !• 

So,  I  shall  replace  cosy  by  1,  assuming  that  Yl,  is  quite  small. 

I  assume  intuitively  that,  in  the  zone  Z0,  9  =  x  and  9’  =  ^ ;  in 
particular,  that  9£  = 

If  the  variation  of  t'  is  very  feeble  in  this  zone  —  which 
occurs,  in  general  —  we  can  substitute  for  9,,*which  gives, 
with  a  great  approximation,  in  the  zone  Z0, 

p  =  6*  +  t£. 
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Furthermore,  by  substituting  the  approximate  expression  A(t) 
for  p,  we  obtain  formula  (70).  This  formula  puts  in  relief  the 
points  of  inflection  of  the  curve  6(t)  —  an  important  fact. 

Then,  I  again  take  the  approximate  expression  (69)  for  6. 

This  approximation  will  be  good  if  |  p/x^  J  remains  small  in  the 
zone  Z0.  In  particular,  if  p0  /  0,  the  ratio  |pe/  x *,  j  will  have 
to  be  small  for  the  approximation  to  be  valid. 

In  general,  the  smaller  |  p/  x  ^ |  is,  the  better  the  approxi¬ 
mation  will  be. 

In  these  questions,  we  must  always  try  to  see  the  conditions 
which  concern  a  very  good  approximation.  If  great  difficulties 
sometimes  keep  us  at  a  little  distance  from  the  conditions  re¬ 
quired  for  a  good  approximation,  we  shall  then  knov/  that  the  method 
is  less  sure,  end  shall  have  to  introduce  some  prudence  in  the 
terms  of  the  conclusion. 

For  example,  when  I  say  that Q  or  w  is  a  large  number  with 
respect  to  certain  parameters,  this  fact  must  be  discussed  and 
verified.  A  complete  discussion  of  the  approximation  would  involve 
a  classification  of  projectiles  according  to  their  diameter,  their 
ballistic  coefficient,  etc.,  and  it  wrould  be  seen  that  such  an 
approximation,  which  is  excellent  for  certain  calibers,  is  weak  for 
others.  I  content  myself  here  with  quite  general  indications. 

In  my  study  of  the  system  (35),  I  always  supposed  P  constant, 
and  substituted  P0  for  P.  That  requires  some  enlightenment. 

According  to  the  preceding, 

(P’/Po)  =  -  n/10,'  (71) 


where  n  is  quite  close  to  1,  but  I  am  not  more  precise. 

Now,  suppose  V0  comprised  between  500  and  1500  m/s.  To  give 
an  indication  that  is  probable  but  not  precise,  I  shall  say  that  Q 
could  be  comprised  between  1000  and  3000.  Moreover,  taking  prob¬ 
able  values  for  the  ratio  A  /B  ,  we  could  say  that  co  is  comprised 
among  such  numbers  as  70#,  100#,  and  150#. 

In  these  conditions,  the  zone  Z0  would  correspond  to  a  dura¬ 
tion  of  2/70,  2/100,  qnd  2/150  sec. 

Now,  let  us  study  the  ratio  b1 /b  in  the  zone  Z0. 


Let  us  first  suppose  60  =  0.  We  have: 
8*  *  b'0,  b  *  8»t, 

6i/&  *  1/t. 


(72) 
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In  the  zone  Z0,  this  derivative  is  extremely  large 


Now,  let  us  suppose  8o>0,  and  extremely  small.  This  concerns 
almost  perfect  departure,  and  I  suppose  essentially  that  the  ratio 
8»/60  s  |  [  /60  is  a  large  number. 

In  the  zone  Z0, 

6t  «  6^,  8  *  60  +  6»t. 

Set 

’  8;t.=  L80. 

For  t  =  0,  * 

L  =  0.  • 

At  the  end  of  the  zone  Z0, 

t  =  2/70  or  2/100  or  2/150, 
approximately,  and,  e.g.. 


L  =  2,  3 ,  4, 

Write 

*  5° 

8  60(1+L) 


(73) 


Evidently,  this  ratio  is  quite  a  large  number,  a  little  less 
than  5*/80,  in  the  zone  Z0,  if  co  is  a  large  number. 

In  all  cases,  the  ratio  6’/&  is  much  larger  thanjPVPjin 
the  zone  ZOJ  if  co  is  a  large  number. 

In  the  zone  Ze,  the  product  P6  is  increasing,  altho  P  may 
be  decreasing;  the  variations  of  P  are  insignificant  with  respect 
to  those  of  6,  in  general. 

To  avoid  all  equivocation,  I  must  repeat  that  I  do  not  examine 
all  possible  cases,  when  I  refer  to  the  numerical  value  of  2  ,  o>, 

or  L.  That  would  require  a  much  more  detailed  discussion. 

12.  SUMMARY  STUDY  OF  AN  ALMOST  PERFECT  DEPARTURE.  I  have 
defined  the  almost  perfect  departure.  In  particular,  (M/6) 0  =  0, 
so  that,  in  the  zone  Z0,  we  shall  put  the  criterion  (30)  in  the 
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form: 


s  +  >2  . >  0,  (74) 

or  else: 

S  +  4^2  >  0.  (75) 

When  pr  is  negative,  we  must  seek  the  maximum  of 
Jp’/M  and  the  criterion  takes  the  form: 

8>  45 2  f -5— 1  max»  (76) 

I  recall  that 

S  =  (  A  :>i)2  -4 B  Rki_. 

Furthermore, 

if1!  max  ««[(§)  max  *P«  )  ,  (77) 

so  that  we  ought  to  seek  the  maximum  of  q/6. 

Set 

<§>aax  =  <8  +  1)  p». 

so  that 

(79) 

the  form: 

(80) 

curve  of  q,  or  B(t),  would  have 
to  be  made  in  the  case  of  almost  perfect  departure.  To  abridge 
the  discussion,  which  would  require  quite  long  developments,  I 
shall  only  make  a  sketch,  by  examining  two  particular  cases. 

First  case.  I  suppose: 


if  1 


max*  cosP0, 
and  the  criterion  will  take 
S  >  4  5s(Rkf)0, 

for 


coP  = 


Rkt 

B 


A  general  study  of  the 
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Then: 


Qq  s  r«  3  0,  Ho  =  0. 


Moreover,  the  ratio  |  p0Aj»[  Is  very  small,  which  will 
permit  me  to  keep  the  method  of  approximation  previously 
adopted.  Here,  as  In  the  case  of  perfect  departure, 

11m  M  n 
t=0  6  = 

but  with  pQ  /  0;  so,  this  is  an  almost  perfect  departure. 

Calculate  the  solutions  p  and  q  of  the  system  (35) , 
which  vanish  at  the  origin;  they  are  represented  by  the  curves 
already  constructed  (Fig.  7).  To  pass  from  p  to  A(t)  and 
from  q  to  B(t),  it  is  sufficient  to  add  respectively  the  periodic 
functions  Wi  and  W2 .  These  functions,  defined  by  formulas  (60) 
and  (61),  have  here  the  simpler  form: 

Ui  =  \,  cos  cot,  (81) 

U2  =  sin  cot.  (82) 

Consequently,  the  curve  of  q  will  be  raised  in  the  interval 
from  T  =  u  to  T  =  2w,  and  the  curve  of  p  will  be  raised  in  the 
interval  from  T  =  */2  to  T  =  3*/2. 

The  solid  line  (Fig.  12)  represents  the  curve  of  q,  and  the 
dotted  line  represents  the  curve  B(t) .  V/e  ought  to  draw  the  tan¬ 
gent  0E  to  the  latter  curve;  it  is  evident  that  the  tangent  is 
raised.  At  the  same  time,  the  straight  line  0e»,  which  represents 
the  first  approximation  of  6,  will  be  lowered  a  little.  So,  the 
maximum  of  q/6  is  augmented  by  the  presence  of  the  non-vanishing 
parameter  p0.  .  - 

Therefore,  if  we  take  M.  de  Sparre’s  coefficient  s  -  1/3 
in  the  case  of  perfect  departure ,  it  will  be  possible  to  have 
s  >1/3  in  the  case  of  almost  perfect  departure. 

Second  case.  Take: 

V0  =  1200,  g  cos  a  =  l,a  *  84 

Then 


IU  I  =  1/1200 


We  shall  take 


80  =  1/120,000. 

which  gives 

50/  |  t»  |  =  1/100. 

The  first  approximation  of  6  will  be 

6  =  1200  (l5o  +  ^  * 


Altho  a)  is  unknown  a  priori ,  we  know  that  its  numerical 
value  will  be  quite  large  if  the  diameter  of  the  projectile 
•is  not  very  large.  So,  suppose  the  caliber  medium,  e.g.: 

80  mm,  100  mm,  140  mm,  ....  Represent  the  value  of  6  for 
T  =  «  by  5  ,  and  the  value  of  8  for  T  —  2*,  which  corresponds 
to  the  endmof  the  zone  Z0,  by  8  .  For  example,  if  jv/co  were 
near  1/200,  we  would  have: 


8  * 
m 


35 

2 


n  I 


fcf 


n 


28, 


If  tt/ a)  were  near  1/100,  we  would  have: 


280,  8 


n 


36 


o  ♦ 


Since  the  variation  of  6  is  not  large  in  the  zone  ZQ,  I 
shall  be  able  to  make  an  approximation  which  is  somewhat 
rougher,  but  rapid  and  convenient.  I  shall  take  a  mean  value 
of  P6  in  Z0 »  the  value  at  the  instant  T  =  #.  The  fact  that 
co  is  not  knov/n  is  unimportant,  as  we  shall  see.  So,  the  mean 
value  of  P5  will  be  Po^,  very  approximately. 

So,  I  integrate  the  system  (35) ,  by  setting: 
fx  =  0)Po6m,  fa  =  to  ~  Po6m* 


We  know  the  solutions  p  and  q,  which  vanish  at  the  origin 
and  we  know  that; 

f  p  *max  *  ? qmax  *  2P°6m  * : 

To  obtain  the  maximum  of  q/6,  I  ought  to  draw  a  tangent 
to  the  curve  thru  the  point  0”.  I  can  assume, 'as  a  new  approx 
imation,  that  the  point  of  contact  E  is  confused  with  the 
summit  2  of  the  curve.  By  the  hypothesis  which  was  made  about 
the  numerical  value  of  A  /p,  the  maximum  of  q  occurs  for 
T  *  *. 
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So: 

®max*2p<>-  (8*) 

This  approximate  result  is  independent  of  the  numerical 
value  of  8  ,  and  that  is  what  makes  the  approximation  easy. 

In  the  actual  case,  with  this  approximation,  s  *  1. 

How,  suppose  the  departure  almost  perfect,  i.e.,  p0, 
qQ,  r0  non-vanishing  and  small  enough  for  the  adopted  approx¬ 
imation  to  be  valid. 

From  the  point  of  view  of  the  approximation,  we  must 
examine  the  numerical  value  of  P0.  The  approximation  v/ill  be 
good  if  P0  is  quite  small,  i.e.',  if  the  ballistic  coefficient 
is  very  small.  That  results  from  the  examination  of  the  formulas 
(S3).  ,  . 

Suppose  the  initial  values  p0  and  qp  negative.  The  curve 
of  q  being  represented  by  the  solid  lire  (Fig.  13),  the  curve 
B(t)  will  be  represented  by  the  dotted  line;  the  curve  will  be 
raised  in  the  region  of  the  summit  2  ,  and  the  maximum  of  B  will 
be  able  to  have  the  same  abscissa  as  the  maximum  of  q.  Then: 

<5>max“(2  +  h>P«>  h  >0>' 
i.e. ,  s  >  1. 

Since  this  method  of  approximation  is  not  very  precise, 
the  value  of  s  is  not  very  well  determined;  but  I  have  shown 
intuitively  that  a  minimum  initial  perturbation  can  make  the 
numerical  value  of  s  greater. 

Since  the  criterion  relative  to  Q  is  represented  by  formula 
(80)  at  the  origin,  we  see  quite  simply  that  minimum  irregulari¬ 
ties  at  the  departure  might  increase  the  coefficient  s. 

Moreover,  if  the  artillery  adopts  a  coefficient  s  greater 
than  M.  de  Sparre's,  that  signifies  that  it  foresees  quite 
strong  initial  perturbations*. 

I  have  examined  one  of  the  numerous  aspects  of  the  problem, 
and  I  shall  complete  these  researches  later.  . 


Translated  by 


\  H.  P.  Hitchcock 

*&,  Sugot,  Theoretic  Exterior  Ballistics,  Mem,  de  l*art.  fr . ,  Vol.  VI, 

1927. 
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